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1. Introduction 

We present a method \^ for improving the estimation of quark propagators between different 
domains of the lattice. Our method turns out to be similar to that of "maximal variance reduction" 
(MVR) However, it contains the advantage that one can work directly with the chosen lattice 
Dirac operator. In the following sections, we present our method and some first results for static- 
hght mesons, where the Chirally Improved (CI) [^], light quark propagator is calculated with the 
improved estimator. 



2. The method 

Decomposing the lattice into two distinct regions, the full Dirac matrix can be written in terms 
of submatrices 

M=f^"^^^V (2.1) 

yM2i M22J 

where Mn and M22 connect sites within a region and Mi 2 and M21 connect sites from the different 
regions. We can also write the propagator in this form: 



M-^=P= I " " " I . (2.2) 
\P2i P22) 

We consider a set of random sources, x" (n = l,...,N), and the corresponding resultant vectors, 
Tj" = Px", to derive useful expressions for our technique. Reconstructing the sources in one region, 
X", from the solution vectors everywhere, T]", we may write 

X",=M7]"=Mn7l'l+Mnri^. (2.3) 
If we now apply the inverse of the matrix within one region, we have 

Mn'Xi =V"+ Mn'MnTl^ . (2.4) 

This can be solved for rj" and substituted back into the naive estimator of the propagator between 
the two regions (repeated source indices, n, are summed over): 

Pn ^ ^v-af 

^ ^[MnHXi-Mi2Tl2)]xf 

« -i(M^/Ml2T]f)z2"^ (2.5) 

where in the last line we eliminate the first term due to the fact that we expect lim^v^oo XiX-i ~ ^^ 
Writing out the full expression, we obtain 

Pn « -^M;^'MnPX2xf 



Pll P\2 



-M^,'MnP22, (2.6) 
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Figure 1: Depiction of one of the boundaries and the surrounding sources (x) which we use to construct 
the two estimates of the CI quark propagator between regions of equal volume. Colors and shapes indicate 
which source region contributes to the signal in the resultant vectors. The y/'s are calculated using only one 
of the source regions, while the 0's use both. 



where the second hne is an exact expression, showing that one can relate elements of different 
regions of P = M^' via the inverse of a submatrix of M. This is nothing new. After all, P22 is the 
Schiir complement of M^^^ . But the lesson learned up to this point is that we need no sources in 
one of the two regions. 



Looking again at Eq. ( |2.6[ ), one can see that we need not make the approximation P22 ~ 
j^PX'iX'i ■ Instead, we can place the approximate Kronecker delta between the M\2 and P22'- 

Pn - -^M^'MuXiXfP 

^-^{Mn'MnX2)iY5Pr5X2f 

«-iv/J'(/.f, (2-V) 

where we have used the 75-hermiticity of the propagator. One can see from the form of the vector 
Yi in the next to last line that we only need sources which "reach" region 1 via one application of 
M. Also, one can use all points in one region for the source and all points in the other region for 
the sink. 

For our first attempt of using this method, we use equal volumes for the two regions and place 
sources next to the boundaries in both regions (see Fig. [I]). Although this choice may not be ideal, 
we do perform inversions for each spin component separately (spin dilution) and, since our sources 
occupy all relevant time slices surrounding the boundaries, we actually obtain two independent 
estimates of the quark propagator between the two regions: 

-^Ym'^Pn « -^750IV/75. (2.8) 
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oper. 




QO{r,Di)q' 


S 


o-,i- 




P_ 


0+,l+ 




P+ 


l+,2+ 






l-,2-,3- 


Qy,{D]-Dl)q' 



Table 1: Static-light meson operators. 



So our method is very similar to that of MVR, except for the fact that we can work directly 
with M, instead of M'M. This is better since it is less problematic to invert M due to it having a 
better condition number than M^M [§]. Also, the sources need only occupy enough time slices to 
connect them with the other region via M, rather than M^'M. These are the same number of sources 
for Wilson-like operators, where M^M, Uke M, only extends one time slice. However, for many 
other improved operators (like CI) this can reduce the number of necessary source time slices by a 
factor of 2. On top of this, although the nature of the lattice Dirac operator may dictate what is the 
ideal domain decomposition, it does not otherwise restrict the choice of regions or the use of this 
method (e.g., it is even possible to use the Overlap operator). 

For expressions and first results relevant to propagators which return to the same region (e.g., 
closed, or loop, propagators) we point the reader to our lengthier publication on the subject |[T]] and 
move on to our application for the open propagators. 

3. Static-light mesons 

For our meson source and sink operators, we use bilinears of the form: 

QO(r,Di,D\S)q, (3.1) 

where 5 is a gauge-covariant (Jacobi) smearing operator [Q] and D is the covariant derivative. For 
our basis of light-quark spatial wavef unctions, we use three different amounts of smearing and 
apply 0, 1, and 2 covariant Laplacians to these: 

q' = Ssq , D'-Snq , D^D^Sieq , (3.2) 

where the subscript on the smearing operator denotes the number of smearing steps; all are ap- 
plied with the same weighting factor of Kg^ = 0.2. So we have a relatively narrow, approximately 
Gaussian distribution, along with wider versions which exhibit one and two radial nodes, due to the 
application of the Laplacians. We point out that, thus far, we have not altered the quantum numbers 
of the meson source since both the smearing and Laplacians treat all spatial directions the same 
(i.e., they are scalar operations). In order to create mesons of different quantum numbers, we use 
these light-quark distributions together with the operators shown in Table |] (see, e.g., Ref. [§]). 
Inserting the estimated and static propagators in the meson correlators we have 



Qj{t) = {0\{QOjq)t (^0;G)o|0> 



1 + 74'"^ 



} . (3.3) 

{u} 
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Figure 2: Effective masses for the static-light mesons on the quenched configurations, am^ — 0.08, a^' « 
1330 MeV, L « 1 .8 fm. The horizontal lines represent M ± Gm fit values for the corresponding time ranges. 
Dashed lines indicate fits for which we adjust the minimum time for systematic error estimates. 

The static quark is propagated through products of links in the time direction and has a fixed spin 
(1 + 74)/2. The estimated propagator P^^^f^^ is of the form of Eq. (2^). Thus, all points within 



region 1 {N^Nt/2 of them) can act as the source location x, just so long as t is large enough to 
have the sink location x + t4m region 2. Note that we now have subscripts on the source and sink 
operators to denote which light-quark distribution is being used. We create all such combinations 
and thus have a 3 x 3 matrix of correlators for each of the operators in Table |^. 

Following the work of Michael and Liischer and Wolff [^, we use this cross-correlator 
matrix in a variational approach to separate the different mass eigenstates. We must therefore solve 
the generalized eigenvalue problem 

C(Ov«=A«(?,?o)C(?o)v«, (3.4) 
in order to obtain the following eigenvalues: 

A«(f,fo) - e-'^'[l + 0{e-'^')] , (3.5) 

where is the mass of the ^th state and AM^ is the mass-difference to the next state. For large 
enough values of t, each eigenvalue should then correspond to a unique mass state, requiring only 
a single-exponential fit. 

Variational approaches have seen much use recently in lattice QCD, especially for extracting 
excited hadron masses and we point the reader to the relevant literature in [^. 

We create our cross-correlator matrices on two sets of gluonic configurations: 100 quenched 
and 74 dynamical, each with 12^ x 24 lattices sites. The quenched configurations have a lattice 
spacing of a « 0.15 fm (a^^ 1330 MeV) and a spatial extent of L ?a 1.8 fm. The dynamical set 
[|] has 2 flavors of CI sea quarks (with M^^^sea ~ 500 MeV), a « 0.115 fm (a^^ 1710 MeV), 
and L PS 1.4 fm. We use 12 random spin-color vectors as sources for the Ught-quark propagator 
estimation. Spin-diluted, this gives us 48 separate sources for the inversions (one in the full volume. 



<p, and two in the subregions, see Eqs. ( [2.7| ) and ([2.8|)). We perform inversions for 4 different 
quark masses: am^ = 0.02, 0.04, 0.08, 0.10. 

After extracting the eigenvalues, we check for single mass states by creating effective masses. 
A representative sample of these, along with single-elimination jackknife errors, are plotted against 
time in Fig. ^ In each case one finds values from the first three eigenvalues. The horizontal lines 
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Figure 3: Physical mass splittings (M — Mis) a function of the quark mass for the quenched (left) and 
dynamical (right) lattices. The vertical lines denote the chiral limit (nig 0) and the strange quark mass 
(mj). The left-most values are the (linear) chiral extrapolations. 



signify the M ± Om values which result from correlated fits over the corresponding range in time. 
We require that at least three effective mass points display a plateau (within errors) and that the 
eigenvectors remain constant (again, within errors) over the same range before we perform said 
fits. 

Performing fits for all quark masses, we next take a look at the mass splittings (M — Mis) a 
function of the quark mass. These are plotted in Fig. ^, along with the chirally extrapolated results 
(niq 0). We use simple linear fits to perform these extrapolations. 

In Table § we present the results for the chirally extrapolated {B mesons) and strange-quark- 
mass interpolated (B^ mesons) mass splittings. We include the statistical errors from the fits in the 
first set of parentheses. For fits where the effective-mass plateaus are not immediately clear (e.g., 
fits represented with dashed lines in Fig. ^), we move the minimum time of the fit range out by 
one to two time slices and observe the subsequent changes in M it Gm, as compared to the previous 
values. The differences from the old values aie reported as systematic errors; these appear in the 
second set of parentheses. For a discussion of these results, we refer the reader to our lengthier 
report 

One thing is clear though: due to the improvement of the light-quark propagator estimation, 
and our subsequent ability to use half the points of the lattice as source locations, we have greatly 
improved our chances of isolating excited heavy-Ught states. In an earlier study of heavy-light 
mesons using wall sources on the same quenched configurations, we were barely able to see the 
25 state, let alone the excited states in any other operator channel. Also, there we used NRQCD 
for the heavy quark; this should only boost the signals since the heavy quark can then "explore" 
more of the lattice through its kinetic term. It is obvious, however, that we have much better signals 
now since we are able to see excited states in every channel {2S, 35, 2P^, 2f+, and 2D±) on the 
quenched lattice. 
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state 


7^ 


M-Mi5(MeV)forS/S^ 






Nf = 0,L^l.Sfm 


A^/ = 2,L?a 1.4fm 


25 


(0,1)- 


712(14)/ 675(10) 


717(69)7 665(45) 


35 


(0,1)- 


1265(40)(+?3o) / 1220(30)(+20) 


1640(44)(+5f,„) / 1560(45)(^35q) 


1P_ 


(0,1)+ 


350(46) / 384(20) 


243(51)7 330(34) 


2P_ 


(0,1)+ 


971(49)(;^°)/923(30)(+^°) 


-/- 


1P+ 


(1,2)+ 


446(15) 7 424(10) 


341(82)7 363(55) 


2P+ 


(1,2)+ 


1028(28)(+i^'') / 993(20)(+i^'') 


1001(80)(+i30)7 930(75)(iOo) 


1D± 


(1,2,3)- 


808(27)(;i,;)/773(17)(+^„) 


-7- 


2D± 


(1,2,3)- 


1183(97)(;}3«)/ 1188(68)(;™) 


-7- 



Table 2: Mass splittings for our B I mesons. The first number in parentheses is the statistical error The 
second set (if present) are the additional systematic errors which result from adjustments to the minimum 
time of the fit. 
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